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Sequential single-electron charging is observed in InAs Self- Assembled Quantum Dots using ca- 
pacitance spectroscopy. In this system, the Coulomb energy is smaller than the inter-level energy 
spacings due to the quantum confinement and both effects can be separately identified. A theoret- 
ical model is proposed for this system and the capacitance experiments were devised in order to 
experimentally observe the effects of Coulomb interaction between electrons on the dots. The effects 
of inter- and intra-dot Coulomb interaction have been observed in the capacitance spectra. A good 
agreement between the proposed model and experiment is achieved. 



I. INTRODUCTION 

Quantum effects and Coulomb effects in zero dimen- 
sional systems have been extensively studied in the past 
years due mainly to thej-pew and interesting phenomena 
these structures exhibit^. The degree of importance of 
each of these effects will depend on the size of the ob- 
jects that form the system, since Coulomb effects should 
diminish linearly with the inverse of the dimensions of 
the system whereas quantum effects diminish quadrat- 
ically. Artificially produced systems, such as those de- 
fined by electron beam lithographya and other pattern- 
ing techniques as well as in structures with split gate 
geometryQ, are usually constrained by technological lim- 
itations. The ultimate size one obtains with such tech- 
niques is still too large to make quantum energies compa- 
rable to thermal energies beyond the mK range, and as 
a rule Coulomb energies prevail in these systems. Nev- 
ertheless, systems where the growth kinetics naturally 
produce small structures offer the possibility of decreas- 
ing even further the dimensions. As an example, is- 
landing jin highly strained heteroepitaxial coherent III-V 
systemso produces structures with dimensions were quan- 
tum effects are routinely .observableoQ, and furthermore, 
at elevated temperaturesLM In this work, the InAs self- 
assembled quantum dots grown on GaAsLrEl are inves- 
tigated with capacitance spectroscopy. Previous works 
on this system using the capacitance technique were fo- 
cused onrthe magneto-optical propertied 9 ] and electronic 
structuret2l of the dots. The focus on the present work is 
on the electrostatic interactions in this system, where we 
analyze the effect of having a finite number of electrons 
in a single dot and the adjoining consequences on the 



capacitance spectra of having their correspondent image 
charges in a nearby electrode and electrons in neighbor- 
ing dots. It is demonstrated that the static character- 
istics show the independent loading of single electrons 
in arrays of about 10 9 dots. We successfully predicted 
and simulated the effects of the Coulomb interaction be- 
tween electrons in the dots with the gate and electrons in 
other dots. The structure of this paper is as follows: in 
the experimental section the growth and the capacitance 
experiment will be described, along with the identifica- 
tion of the components inherent to each capacitance spec- 
trum. The theoretical model section describes the model 
utilized to predict and simulate the capacitance spectra 
and the electrostatic effects, and the results and analysis 
section presents the measured data and the simulation 
results. 



II. EXPERIMENT 

The samples studied in this work were grown on a 
Varian GEN-II MBE system on semi-insulating [100] ori- 
ented 2" GaAs substrates. The conduction band diagram 
of a sample is shown schematically in Fig. |l|. Tho^ewth 
procedure was the same as used in other worksffllI2E3, 
consisting on the oxide desorption under 10 -5 Torr AS4 
flux and subsequent growth of a 1 micron undoped GaAs 
buffer layer at 1 ML/s at 620 °C. We also used 40 periods 
of a 2 nm x 2 nm AlAs/GaAs short period superlattice 
(SPS) to trap defects and smooth the surface prior to 
the deposition of a 20-100 nm GaAs undoped layer, 20- 
80 nm thick 2 x 10 18 cm~ 3 n-doped back contact and a 
undoped GaAs tunneling barrier with thickness d. The 
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FIG. 1. Sketch of the band diagram of the quantum dot 
array sample in the direction perpendicular to the dot and 
gate planes. 



growth temperature is now reduce to 530 °CE3 and InAs 
is deposited, with the substrate fixed with the llO direc- 
tion aligned with the axis of the In cell. By doing that, 
we could produce a variation on the deposited thickness 
of InAs of up to 10%, assuming an ideal point source for 
the effusion cell. This procedure allowed us to obtain a 
the dot density ranging from in one side of the wafer 
up to 2 • 10 10 cm" 2 , as estimated by TEM and capaci- 
tance measurements. The growth was then resumed by 
depositing GaAs at 530 °C with the substrate now ro- 
tating. After 5 nm of GaAs thickness, the temperature 
was increased to 600 °C remaining at this value for the 
remaining part of the growth. A 25 nm thick GaAs was 
the deposited, and another AlAs/GaAs 3 nmx 1 nm SPS 
was grown before capping the sample with a 5 nm thick 
GaAs layer. 

When analyzing the experimental data, one should 
take into account that the potential in the dot plane V 
is not equal to the voltage on the gate V g . If the dots 
are located at a distance d from the back gate and the 
distance between the front and back gates is t than 
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We call the ratio L a lever-arm coefficient, it varies from 
6 to 7 in different samples. We kept values of L similar 
for all samples with different d by varying the thickness 
of the top AlAs/GaAs SPS. 

Samples with tunneling barriers thinner than 20 nm 
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FIG. 2. Dot capacitance vs gate voltage for the samples 
with the distance d — 450A between dot and gate planes. 
Circles show the experimental data, solid line gives the fit by 
the sum of the linear background and two Gaussian peaks 
(shown by dashed lines). 



were not used since the strain field produced by the 
dots extends into the substrate for approximately this 
amount, as qualitatively inferred by changes in contrast 
in cross-section TEM micrographs. The thickness of the 
back-contact had also to be increased so that the strain 
field would not deplete the regions beneath the dots and 
therefore producing effectively different d values. 

Metallic disks were evaporated after contacting the 
back contact layer using standard procedures for making 
ohmic contacts in n doped GaAs. The photolithographi- 
cally defined disks consisted of a thin layer of Cr (5 - 10 
nm) followed by a 200 nm thick Au layer, with diameter 
from 150 /ini to 350 /im. The diodes were bonded and 
their capacitance was measured at 4.2 K in a liquid He 
immersion cryostat. To obtain the capacitance voltage 
characteristics of the devices, we used a lock-in amplifier 
at frequencies varying from 40 Hz to 30 kHz. The ac 
amplitude was kept smaller than 1 meV, after conversion 
into energy by the lever arm coefficient. Depending on 
the thickness of the tunneling barrier d we used higher or 
lower frequencies, choosing the appropriate value to have 
the loading and unloading of the dots within one period 
of the ac frequency, i.e., / < 2ttRC, where R is the tun- 
neling resistance and C is the capacitance determined by 
the area of each dot and the distance d from the back 
contact. 

In Fig. ^ we show the measured capacitance-voltage 
characteristics of the sample A with d = 450A. One can 
see the two peaks, corresponding to the loading of the 
ground and first excited states of the dots. The first peak 
is split into two; this splitting results from the Coulomb 
interaction of the electron in the same dot. In this paper 
we analyze only this split peak, related to the ground 
state in the dots. 
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The following problem complicates a direct comparison 
of the experimental results with theoretical capacitances: 
the experimentally measured capacitance C OX p contains 
a large background capacitance Cb, to which the capaci- 
tance of the dots Cdot is a small correction: 
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a 



dot- 



We resolve this problem by assuming a linear background 
capacitance Cb and extrapolating it using the low- voltage 
part of the CV curve. For the sample A we have also mea- 
sured the capacitance at high frequency of the ac voltage. 
In this case the electrons do not have the time to tunnel 
to and from the dots during one period of the ac voltage, 
and the dots essentially do not contribute to the capaci- 
tance, whereas the background capacitance has no intrin- 
sic slow processes, and does not depend on the frequency 
in a wide frequency range. Thus measured background 
capacitance shows a linear voltage dependence for the 
voltages where at low frequencies the dot capacitance is 
observed. We have also used a background capacitance 
in the form of the capacitance of a depletion layer, which 
is formed in the doped back contact: 
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Here the indices i, j number the quantum dots, is 
the occupation number of the dot i; in this paper we 
only consider the ground state on the dot, so rii can be 
0, 1, and 2. The first term represents the non-Coulomb 
disorder due to fluctuations of size quantization energies 
ipi in the dots. We assume here a Gaussian distribution 
for ip: 
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The width a v describes the magnitude of disorder. It 
is a fitting parameter in our model. The second term 
in the Hamiltonian Eq. is the charging energy of the 
dot (this term is responsible for Coulomb blockade-like 
effects). Ci is the capacitance of the dot, defined so that 
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where Vd = e 2 Nd/K, Na is the density of the positive 
charge in the depletion layer, n is the dielectric constant, 
Vo is the voltage of flat bands, and S is the area of the 
sample. In the range of V where the capacitance peak is 
observed (and for smaller voltages), this function can be 
considered as a linear one. All three ways of extracting 
the background capacitance lead to practically the same 
results. 

After the background is subtracted, the capacitance 
correction due to quantum dots Cdot can be very well 
described by a sum of two Gaussian peaks with centers 
at voltages V\ and Vi and widths a\ and &i , respectively: 
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This fit reflects the fact that the areas of both peaks Qq 
should be equal, since they give the total charge of all 
dots with one electron per dot. The fitting was done by 
weighted explicit orthogonal distance regression using the 
software package ODRPACKcJ. The resulting decompo- 
sition of the experimental curve into linear background 
and two Gaussian peaks is shown in Fig. @. 



III. THEORETICAL MODEL 

We now present a theoretical description of the quan- 
tum dot array in the framework of the following Hamil- 
tonian: 



kC 



is the energy of Coulomb interaction of two electrons in 
the dotfij. The product m(rii — 1) takes into account 
that electron does not interact with itself. In our model, 
we take Ci — C to be the same for all dots, and use it 
as another fitting parameterli-3. The third term in the 
Hamiltonian is the attraction energy between the charge 
of the dot and its image in the metallic electrode at a dis- 
tance d from the dot plane. We have considered only the 
attraction to the closest metallic electrode, which is the 
back gate in our samples. The front gate is much further 
away from the dots, and has no significant effect on the 
electrostatics near the dots. The last term introduces the 
effect of the gate bias V (scaled by the lever-arm ratio 
Eq. (H). 

The fourth term of the Hamiltonian Eq. (g) takes into 
account the interaction between the charges on different 
dots. The interaction potential V(r) includes the screen- 
ing of the Coulomb interaction by the gate and has the 
form: 
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Here we again neglect the effects of the second, remote 
gate. We should also mention that we have assumed 
the point charge Coulomb interaction of the dots with 
the gate and with each other The latter assumption is 
well justified when the dot size is much smaller than the 
distance between the dots, which is the case in all of 
our samples. For the dot-gate interaction, a trivial cal- 
culation shows that the interaction of an electron wave 
function in a disk of size a with the gate at a distance d 
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differs negligibly from that of a point charge when d> a. 
This condition also holds for all of samples for which we 
present the results here. 

The thermodynamic properties of the many-electron 
system with the Hamiltonian Eq. (|^) in the regime when 
the Coulomb interactions are essential are rather com- 
plex, and a comprehensive solution is possible only by 
means of numerical modeling. The modeling consists of 
the following steps: 

1) Generation of dot positions and energies. The given 
number of dots N are placed inside a square of the 
size L. The dot density N^ot = N/L 2 for our sam- 
ples is about 10 10 cm~ 2 , and we can simulate systems 
of N = 2500 - 10000 dots. The positions of the dots 
are random, with one restriction: no two dots can be 
closer to each other than the minimum dot separation 
r m ; n . The results weakly depend on r m ; n ; we can deduce 
its approximate value r m i n = 250A from the TEM data 
for our samples. Also at this step we generate Gaussian 
random energies ifi. 

2) Monte-Carlo step. On this step we use the Metropolis- 
type Monte-Carlo algorithm to calculate an average den- 
sity of electrons on the dots n — (rii) for a given gate 
voltage V g and temperature T. In this paper we have not 
studied the temperature dependence of the capacitance, 
and the temperature was maintained at T = 4.2 K. We 
employ the parallel version of Metropolis algorithm for 
Coulomb systemslij. The gate voltage is being slowly 
varied across the desired interval, so that the system has 
enough time to reach equilibrium for each value of V g . 

3) Calculations of capacitance. In order to find the ca- 
pacitance Cdot of the dots as a function of gate voltage, 
we compute the derivative It is easy to show that 
the correction to the capacitance per unit area can be 
calculated from the equation 

-£?<'-?>■ W 

The voltage V is connected to the gate voltage V g by 
Eq. (1). This capacitance has a lot of "noise" unless 
the temperature is very high or the number of dots N is 
large, and we get a smooth curve by convoluting it with 
a Gaussian filter with small width. The results of these 
computations are then compared with the experiment. 

The computations were done on the Intel Paragon us- 
ing 128 nodes. It takes about 6 hours to compute one 
C-V curve. 



IV. RESULTS AND ANALYSIS 

We begin by comparing the experimental CV-curves 
with the results of our modeling (Fig. |j). First of all, 
we subtract the background capacitance from the exper- 
imental data as described above. Since the areas of the 
real and simulated devices are different, so the two ca- 
pacitances cannot be directly compared. We therefore 
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FIG. 3. Dot capacitance vs potential in the dot plane 
for two samples with different distances d between dot and 
gate planes: a) - A, d = 450A, b) - B, d = 200A. Cir- 
cles show the experimental data with background capacitance 
subtracted, solid lines give the results of our modeling with 
non-Coulomb disorder a = 8.5meV and dot charging energy 
W = e 2 /nC = 23meV. 

normalize the calculated capacitance so that the areas 
under the curve, which give the total charge accumu- 
lated on the dots at large voltages, are the same. We 
have checked that the area of the experimental capaci- 
tance peak corresponds to two electrons per dot for our 
sample sizes and dot density N^ot ~ 10 10 cm~ 2 , which is 
in good agreement with the TEM measurements of N^ ot . 
We also convert the gate voltages to the potentials in the 
dot plane using the lever-arm coefficient Eq. (|l|) . 

The Fig. (||) shows the measured (circles) and com- 
puted (lines) dot capacitances for two samples with dif- 
ferent dot to gate distances d = 450A and d = 200A. The 
parameters of the modeling were chosen to get a best fit 
for the sample A with d = 450A (Fig. (||a): dot charging 
energy W — e 2 / kC = 23meV and non-Coulomb disor- 
der magnitude a v — 8.5meV. We then computed the 
capacitance for the second sample, B with d — 200A, 
using the same parameters. One can see from Fig. Qb 
that the theoretical and experimental capacitances agree 
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quite well. Since both samples were grown under similar 
conditions and differ only in the thickness of the tunnel- 
ing layers, it is reasonable to assume that the sizes of the 
dots and their fluctuations should be the same for both 
samples. Therefore, the ability of our model to reproduce 
the effect of the change of d on the capacitance, which 
is quite substantial, without any additional parameters 
should be considered as an evidence in favor of the model. 

We now turn to discussing the effects of the Coulomb 
interaction on the dot capacitance. First of all, the 
Coulomb interaction of two electrons in one dot is re- 
sponsible for the structure of the capacitance peaks. In 
the absence of disorder the capacitance would exhibit two 
delta-function peaks separated by the dot charging en- 
ergy. These peaks are located at those values of gate volt- 
age V g i and V g 2 at which first and second electrons, re- 
spectively, enter the dot. The corresponding potentials in 
the dot plane are V\ = V g i/L and V2 — VgijL. The disor- 
der broadens the peaks, however, they are still observable 
in our samples. The charging energy depends first and 
foremost on the dot shape. It has, however, a correction, 
which comes from the interaction of the charges in the 
dot with their electrostatic images in the metallic gate. 
It is easy to show that this correction reduces the voltage 
between the capacitance peaks 6V = V2 — V\ by 
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This correction leads to the observed dependence of the 
peak splitting on the distance d between the dots and the 
gate. 

This picture is further complicated by the interaction 
of electrons on different dots. This interaction gives rise 
to two effects. First, in order to load an electron into 
a dot, one has to overcome the repulsion of all charges 
already on the dots. This effect shifts the position of ca- 
pacitance peaks to higher voltages. The second peak is 
shifted much more than the first one, since the shift is 
proportional to the charge density already on the dots. 
Therefore, the dot-dot interaction increases the splitting 
between capacitance peaks SV. This correction to SV, 
which we call SVmt, is also ^-dependent: the smaller is 
d the weaker is the screened Coulomb interaction and 
the smaller is the correction. One can see that the in- 
dependence of SVi a t and SV gatc are the same - the peak 
splitting increases at large d due to both corrections. A 
very rough estimate of <5Vint can be obtained by calcu- 
lating the average potential created by a random set of 
singly charged dots (with minimal distance between near- 
est dots r min ): 
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The second effect of the interdot interaction also has 
to do with the potential created by the charges on the 
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FIG. 4. Capacitance peaks splitting SV (a) and peak 
widths o 1 and 02 (b) vs the reciprocal dot-gate distance 1/d. 
The experimental values of SV and o\ are shown by circles, 
those of 02 - by triangles. The results of the modeling are 
shown by solid squares. The modeling gives very close values 
of o\ and 02, so only a\ are shown. The dashed line is a guide 
for the eye. The solid line shows the slope of the correction 
SVgate, Eq. (|). 

dots. Since the dots are placed randomly in the plane, 
the potential created by them is also random. For each 
dot, this potential is added to the ground state energy, 
resulting in effectively increasing the total disorder in the 
system. This additional, "Coulomb disorder", is also in- 
dependent. Its order-of-magnitude estimate is provided 
by the dispersion of the screened Coulomb potential of 
randomly placed dots: 



Saf nt = N dot / V(rY d 
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This "Coulomb disorder" will add to the disorder caused 
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by fluctuations of the size quantization energies in the 
dots, resulting in wider capacitance peaks. 

One can easily see the above interaction effects in the 
modeling. For example, the theoretical curve in Fig. ^a 
is obtained for charging energy W = 23mcV and dis- 
order magnitude er = 8.5meV, while the distance be- 
tween the peaks is SV « 24.5meV and peak widths are 
&i ~ 02 ~ 9.5meV. However, experimentally all of the 
effects of the Coulomb interaction are best studied using 
their sensitivity to the screening by the metallic gate. 
All of these effects are qualitatively evident already in 
the Fig. ||, where one can see that for the sample A with 
larger d = 450A the peaks are wider and further apart 
than for the sample B with smaller d = 200A. 

For a quantitative analysis, in Fig. ^ we compare the 
measured and computed values of peak splitting SV and 
widths o i , a% for a set of samples with different dot-gate 
distances d. Again, the dot capacitance C and the non- 
Coulomb disorder a were chosen from best fit for one 
sample with d = 450A, and then the results were com- 
puted for other values of d without any additional ad- 
justable parameters. One can clearly see that the peak 
splitting increases with d, and that this increase is faster 
than the correction SV ga tc (Eq. (0)) alone would give. 
This shows that the dot-dot Coulomb interaction, which 
is responsible for the second correction SV[ n t, is very im- 
portant in the arrays of self- assembled quantum dots. 
The modeling includes both of those corrections, and it 
provides very good description of the experimental data. 
The peak widths also increase with d as a result of the 
Coulomb interaction (see Fig. ^b), again with modeling 
providing reasonably good description of the experiment. 

V. CONCLUSIONS 

We have presented a capacitance spectroscopy study of 
In As self- assembled quantum dots grown on GaAs and 
compared the experimental results with the computer 
simulations. By doing such a comparison for a series 
of devices with different distances between the dot and 
gate planes, we have shown the effect of the Coulomb in- 
teractions, both inside one dot and between the dots, on 
the capacitance spectra. 

We have found out that the Coulomb interaction pro- 
vides two capacitance peaks which correspond to sequen- 
tial loading of two electrons in one dot. The distance 
between the peaks and the width of the peaks contain 
an information about interactions inside one dot and be- 
tween the dots. In our experiments the interdot inter- 
action was weaker than the external (non-Coulomb) dis- 
order. A new, interesting physics could be observed in 
an ordered dot array with weak disorder: then the inter- 
action controls the distribution of electrons on the dots. 
This may change Ihe sign of the derivative as was 
predicted by EfcosO and was observed by Eisenstein, Pf- 
effer, and Westli3 in clean 2D systems. 
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